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ABSTRACT 

We reformulate the heavy quark effective theory in the presence of a residual 
mass term, which has been taken to vanish in previous analyses. While such a 
convention is permitted, the inclusion of a residual mass allows us to resolve a 
potential ambiguity in the choice of the expansion parameter which defines the 
effective theory. We show to subleading order in the mass expansion that phys- 
ical quantities computed in the effective theory do not depend on the expansion 
parameter. 



Submitted to Nuclear Physics B 



* Work supported by the Department of Energy, contracts DE-AC03-76SF00515 and 
DE-FG03-90ER40546. 



1. 1. Introduction 



There has been much recent interest in the limit of QCD in which the mass of 
a heavy quark is taken to be much larger than the characteristic scale Aqcd of the 
strong interactions [1-14]. In this limit it is natural to consider an expansion of the 
QCD lagrangian in inverse powers of the heavy quark mass, as well as kinematics in 
which the heavy quark, even though bound into a hadron, is almost on shell. Such 
an expansion has been used to lowest order in models of heavy hadrons [15,16], in 
QCD sum rule calculations [17-21], and in lattice gauge theory [22-25]. However, 
to put the predictions of this heavy quark effective theory (HQET) on firm footing 
it will be necessary to perform computations to subleading order in the expansion 
in inverse mass 1/mg. But here arises a potential ambiguity, namely in what mass 
should one expand: the pole mass computed in perturbation theory, the MS mass, 
the mass of the physical hadron, or some other parameter? Our intuition tells us 
that nothing physical can depend on this choice, and this is in fact the case. It is 
the purpose of this note to demonstrate in detail how physical matrix elements are 
insensitive to variations in the expansion parameter. 

Throughout this paper we will use dimensional regularization to cut off the 
ultraviolet divergences which arise in the HQET. We will not address the important 
issue of whether power divergences, which appear when a dimensionful regulator is 
employed, require additional nonperturbative subtractions in the effective theory 
[26]. However, the arguments given here are sufficient to show that if such effects 
are important and introduce additional uncertainties into computations, they are 
unrelated to ambiguities in the choice of the expansion parameter which defines 
the effective theory. We feel that it is important to disentangle these two issues. 

The potential ambiguity in the choice of expansion parameter arises because 
in a confining theory such as QCD there is no true "pole" mass which can be 
assigned to a heavy quark. In the absence of such a canonical choice, a variety 
of perturbative prescriptions (e.g., the pole mass to a given order, or the mass 
renormalized at some scale) compete with "physical" prescriptions such as the 
mass of the lightest hadron containing the heavy quark, or this mass minus some 
fixed number. 

In this paper we will show that these various choices correspond in the effective 
theory to various values of a new parameter, the "residual mass" 5m. Section 2 is 
devoted to the construction of this more general form of the HQET. In Section 3 
we discuss, to subleading order in 1/mg, hadronic matrix elements of currents 
containing two heavy quarks. Currents with one heavy and one light quark are 
addressed in Section 4. Section 5 contains a brief summary of our results. 
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2. 2. HQET with a Residual Mass Term 

Let us fix any reasonable choice for the heavy quark mass and denote it by 
rriQ. Each value of mq yields a different effective theory; it is the parameter 
which defines, once and for all, the relationship between the HQET and full QCD. 
Following refs. [9-11], we derive the effective theory as an expansion in I/tuq. We 
rewrite the heavy quark momentum P^ as 

P/j, = rriQV^ + k^ , (1) 

where is the four-velocity of the hadron containing the heavy quark; in the 
mq — > oo limit it is also the velocity of the heavy quark itself. Note that does 
not scale with mq; in fact, it is this property of separating out the part of P^ 
that scales with the heavy quark mass which characterizes a "reasonable" choice 
of rriQ. We also define effective fields hq(x) in terms of the original fields Q(x) by 

hq(x) = exp(img ylv ■ x) Q(x) . (2) 

These fields create and annihilate heavy quarks and antiquarks moving at velocity 
v^. We further specialize to heavy quarks by imposing the condition fthq = hq. 
To lowest order in the 1/mq expansion, the effective lagrangian is obtained by 
matching matrix elements in QCD with those in the effective theory. We find 

£ e ff = hq(iv ■ D)Hq - 5m hqhq = hqv^{\D^ - 5mv^)hQ , (3) 

where = — igA^ is the gauge-covariant derivative. 

The term proportional to 5m corresponds to a residual mass for the heavy quark 
in the effective theory, of a size of order Aqce). In general, such a term arises due 
to an incomplete cancelation of the full theory mass by the field redefinition (2). 
Changes in mq will induce changes in this uncompensated mass: if mq mQ + a, 
then 5m — > 5m — a. Similarly the residual momentum changes according to k^ — > 
k^ — av^. We may write these relations in differential form as 

9Sm = _ 1 9k t _ = _ v ^ (4) 
dm Q ' dmQ M 

The residual mass term has been neglected in all previous derivations based on 
the HQET. This is because mg always has been fixed implicitly, by imposing post 



* Note that the MS mass mQ^rriQ) is not a "reasonable" choice, since in this case contains 
a piece proportional to a s (mQ)mQV^. 
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facto the condition 5m = 0. Clearly there is some choice of mg which accomplishes 
this, even if we cannot compute it. In other words, these derivations have fixed the 
equations of motion of the heavy quark in the HQET to be (if • D)Iiq = 0, rather 
than the more general condition {iv • D — 5m)hQ = 0. Of course setting 5m = 
is as valid a choice of mQ as any other; in fact, for perturbative calculations it is 
by far the most convenient one. Yet nowhere it is forced upon us. Finally we note 
that although it is possible to formulate the HQET with any mass 5m of the order 
of AqcDj if there is more than one heavy quark in the theory it is convenient to 
choose 5m to be the same for all of them, so as not to violate the heavy quark 
symmetry explicitly at leading order. 

Essentially the only modification which we have introduced into the effective 
theory is the change in the equation of motion for the heavy quark, which we may 
write in momentum space as 

v^ikf, - 5m VfI + gAp) h Q = . (5) 

This equation is used in the matching conditions which determine the effective 
lagrangian to sub leading order in 1/mg, as well as in the matching of currents 
in the full theory onto currents in the effective theory. Recall that we may use 
the equations of motion here because the matching conditions are obtained by 
requiring the equality of physical matrix elements in the two theories. 

We have recalculated, to leading logarithmic order in perturbation theory and 
to order 1/rriQ, the operators and currents which arise in the HQET with two 
heavy quarks, which we shall call for concreteness b and c. This corresponds to 
QCD plus a contact term which describes the semileptonic decay of b to c quarks: 

£ fun = b(i IP - m b )b + c(ip- m c )c + cT^bL M + £ light , (6) 

where Aight describes the interactions of the light quarks and gluons. Here is 
the charged lepton current (which includes constants such as Gp and V c b), and we 
have written the left-handed heavy quark current in the general form cT^b. In the 
HQET, £ii g ht is unchanged, while the kinetic and interaction terms are expanded in 
powers of 1/mQ. The strong interaction lagrangian for the heavy quarks becomes 

£kin = h h (iv h ■ D - 5m)h h + -^—h h (iD - 5mv h ) 2 h h 

2m h ^ 

+ a h (n)^-h h a IM ,G^h h + (b->c) , 

where we have ignored an operator whose matrix elements vanish by the equations 
of motion. Such operators do not contribute to physical matrix elements. In lead- 
ing logarithmic approximation, oq(//) = [a s {mQ)/a s (fi)} 9 ^ , where j3 = 33 — 2Nf, 
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and Nf is the number of light quark flavors. The /z-dependence of the coefficients is 
required to cancel the subtraction-point dependence of the renormalized operators 
in the usual way. We emphasize that and m c are fixed (/i-independent) expan- 
sion parameters which define the effective theory. Finally, the effective lagrangian 
is written in terms of renormalized heavy quark fields, which are related to the 
bare fields by a factor Zq 2 = [a s (mQ)/a s (/i)] 4 ^. 

The interaction term in the effective theory is fixed by the matching of the 
currents J = cTb in the full theory onto linear combinations of effective currents. 
Decomposing to order I/tuq the full theory on-shell spinor u in terms of the two- 
component spinor Uh which describes the heavy quark in the effective theory, 



u 



l + ^(^-5m)\u h , (8) 



we obtain the matching conditions at tree level. Additional operators are then 
induced by the renormalization group. To leading logarithmic order, the result is 

Jeff = Co(/i)Qo + Ci(/i)<2i + C 2 (/i)Q2 + C3(/i)<?3 + C^{jj)QA , (9) 

where 

Qo = e^h c Th h , 

Qi = ^ ^—h c {-i1p - Sm)Th h , 
2tti c 

Q 2 = e 1 ^ ^—h c (-\v h ■ *D - 5mv c ■ v h )Th h , 



^—n c {-w h -iJ-omv c -v h )Ln h , ^ 



Q3 = e i< t>-l—h c r(i0-5m)h h , 
2mb 

= e 1 ^ — h c r(iv c ■ D - 5mv c ■ Vh)h h . 
2mb 

We have again neglected operators whose matrix elements vanish by the equations 
of motion. The phase = -(mb^b _m cfc) compensates for the field redefinition 
(2). At tree level, co = c\ =03 = 1, C2 = C4 = 0, while summing the leading 
logarithms yields 



co(^) = Ci(jj) = C 3 (/i) 



CLL 



cm = cm = - l 4 ( r ^- w ) ( ^m) L ln («>w 



(3 \ w 2 - 1 J \ a s (fx) J V «s(i«) J 



n L , : (11) 
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where 

g 

w = v c ■ v h , a L = - [wr(w) - l] , 

1 , (12) 

r(w) = In (w + v w 2 — l) . 

Vffi — 1 

These expressions arise when the transition from QCD to the effective theory is 
implemented at a single scale fh, which in this case is some mass between m c 
and m\y. Improvements such as including the full one loop matching conditions, 
summing the logarithms between m c and to resolve the ambiguity in fh, or 
working to higher order in perturbation theory, are discussed in the literature [27- 
32]. 

From (7) and (10) it is apparent that in the effective theory with a nonvanishing 
Sm the covariant derivative acting on a heavy quark field fiQ(v) always appears in 
the combination 

YD^v) =iDfj,- SmVfj, , iV ^(v) = iD fl + Smv^ . (13) 

Note that the " magnetic interaction" operator in (7) is also of this type, since G ^ v 
is proportional to [Dn,T> u ]. From (4) we may conjecture that matrix elements of 
the operator T> do not depend on the choice of tuq. We shall now demonstrate 
that this is indeed the case. 

3. 3. Hadronic Matrix Elements 

Let us evaluate the hadronic matrix elements of the weak current in the effective 
theory with a residual mass term, repeating the arguments of ref. [13] in this more 
general context. For simplicity, we shall restrict ourselves to the discussion of l/m c 
corrections only. The inclusion of terms proportional to l/m^ is straightforward 
[16]. We will also specialize, for purpose of illustration, to transitions between 
and Z)W mesons. At next-to- leading order in the heavy quark expansion, one has 
to include the 1 / m c corrections to the current and to the hadronic wave functions. 
Matrix elements of the effective current (9) can be evaluated most concisely by 
using a covariant trace formalism [12,13,33,34], 

(D^(v c )\h c Th h | ) = -£oK^)Tr [M D TM B ] , 

( £>« (v c ) \h c [- i^(uc)] Th h | BW (v h ) ) = - Tr [ £> c , v h , M D TM B } 
where pseudoscalar and vector meson states are represented respectively as 

M(v) = - V ^^ 1 5 , M>,e) = v /^F^ 4. (15) 

The universal Isgur-Wise function £o is the only form factor that appears in lead- 
ing order of the 1/mg expansion [4]. At subleading order, there are additional 



(14) 
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functions. The most general Lorentz-invariant decomposition of £^ is 

e(v c , v h , n) = t+(w, /i) (t£ + <) + £_(«;, I*) (vg - <) - ftK AO 7 M • (16) 

The form factors £j, which parametrize matrix elements in the effective theory, are 
independent of the heavy quark masses. They depend, however, on the renormal- 
ization point /i and on the kinematic invariant w = v c ■ v^- From here on we will 
suppress these arguments. 

Performing an integration by parts on the interaction lagrangian and using the 
equations of motion (see ref. [13] for details), we find that in matrix elements 

h c [-w h - i D(v c )]Th h -^(K-5m)(w-l)h c Th h . (17) 

Here A denotes the difference between the mass of the hadron containing the heavy 
quark and rriQ, i.e. in this case A = m^ — m^ = mjj — m c . The relation (17) implies 
that, in matrix elements, we may replace 

Q 2 ^ -L(A-5m)(w-l)h c rh h . (18) 

The matrix elements of Q2 (and Q4, which obeys an analogous relation) are thus 
proportional to the Isgur-Wise function £q- Furthermore, we see that they involve 
the combination (A — 5m), which is invariant under shifts of m c and m^. 

As in ref. [13], the equations of motion and relation (17) may be used to reduce 
the number of independent form factors. We find 

£- = \(K-5m)iQ , 

? (19) 
(w + l)£+ = -(w- 1)(A - 5m) Co - 6 ■ 

The corrections proportional to l/m^, although we do not present them explic- 
itly, involve the same set of universal functions since the second of eqs. (14) can 
be related by an integration by parts to a formula for the matrix elements of 
hcUV^Vb)^. 

The l/m c corrections to the hadronic wave functions come from insertions of 
the subleading operators 



^— h c (YDfh c + acM-p- h c a^G^h c = Ox + a c (fi)0 2 (20) 



into matrix elements of the leading order currents. This gives rise to additional 
unknown functions. In particular, insertions of O2 induce violations of the heavy 
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quark spin symmetry. Extending the trace formalism to these contributions [13], 
we defined Lorentz and CP invariant form factors Xi y i a the time-ordered products^ 



( D«(t; c ) | i / d 4 yT[h c Th h (0), O x (y) ] \ B^(v b ) > 



Tr [M D TM B ] , 



(D^(v c )\i d 4 yT[h c rh h (0),O 2 (y)} \B^\v h ) 



= -±Tr 



(21) 

Hence to order l/m c the hadronic matrix elements of the weak current may be 
expressed in terms of five universal (real) functions of w and the parameter (A — 
5m). As mentioned above, the corrections proportional to l/m^ involve the same 
functions [16]. 

Evaluating the traces in (14) and (21), we obtain the vector and axial vector 
current matrix elements relevant for semileptonic B — > Z)W decays. To order 

l/m c , they are^ 



D\V\B) = v^5^i{coeo« + <) 



+^o(A-5m)(™-l)(< + <) 
4 



.,„. [2 X i-4( W -l)a c X2 + 12a c x 3 ] « + <)} , 

(22) 

(D*(e)\V i \B) = i^m D ,m B e^ a P elv ca v h p \c £o + Co (A - 5m) 

c (23) 

+ ^-6^-*»)(«-i) + ^-[2xi-^»]}, 



f This definition of X2 and %3 differs from that of ref. [13] by the prefactor a c (/z). 
| We correct an error in expressions (3.11) and (3.14) of ref. [13]. 
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( £>*(£) I A" | B ) =^m D ,m B {co «o [e*"(w + 1) - e*-«b <] 

] 



+ 



+ 



+ 



c 2 ^ /-r 



2m c 
2m r 



+ £o (A - 5m) [(w - l)e*" 
e (A-5m)(w-l) [^(w + l) -e*-v h v£] 

+ [2 X i - 4a cX3 ] [e^(w + 1) - } 



(24) 



It is clear that the requirement that these expressions be invariant under redefi- 
nitions of the expansion parameters m c and m^ is equivalent to the requirement 
that the form factors be invariant, i.e. 

d& dxi 







(25) 



dvrtQ dvriQ 

where mq is any of the heavy quark masses. This is not unexpected, since we have 
defined the form factors in terms of matrix elements of operators built from the 
"generalized" covariant derivative T>. It is these form factors, together with the 
invariant combination (A — 5m) , which are observables of the effective theory. 

We recall that by evaluating matrix elements of the charge operator Vq at 
zero recoil one can show that the Isgur-Wise function is normalized at zero recoil, 
£o(l,/-0 = 1, and that two of the subleading form factors vanish at this point, 
Xi(l,/i) = X3(l ) / i ) — 0, if one uses the hadron masses in the prefactors in (22)- 
(24). These conditions are preserved in the effective theory with a residual mass 
term. 

The physical matrix elements must be invariant not only under shifts in the 
expansion parameter toq, but also under changes in the renormalization point \x. 
We can use this fact to deduce the /i-dependence of the universal functions £o> £3 
and Xi from that of the short distance coefficients Cj. To first order in a s we find 
that 



d/j, 

dxi = 
d/i 

d X i 



4a s 
' 3tt 
4a s 
' 3tt 
4as 
' 3tt 



[wr(w) - 1] & 



0,3 , 



[wr(w) - 1] xi - £ (A - 5m) 



r(w) — w 
w + 1 



(26) 



wr(w) + A 



Xi 



2,3. 



These conditions must be obeyed by any model calculation which is sensitive to 
the /^-dependence, such as leading-log improved QCD sum rules. 
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4. 4. Heavy-Light Currents 

The arguments given above apply equally to currents with one heavy and one 
light quark. We have verified that again the only change with respect to the 
case 5m = consists in the replacement of the covariant derivative acting 
on a heavy quark field by from (13). Derivatives acting on the light quark 
field remain unchanged. The study of heavy-light currents is interesting, however, 
since it will allow us to give an operational definition of the invariant combination 
(A — 5m) which does not contain the residual mass explicitly. 

In leading logarithmic approximation, the 1 / mq expansion of the current J = 
qTQ in the effective theory reads 

1 6 

Jes = b in) tf + — £ (/x) QJ , (27) 

where 

Q' = e i( t>'qTh Q , Q , 1 = ^'qnph Q , (28) 

and <p' = —mQV ■ x. The remaining operators in (27) involve the light quark mass 
m q or derivatives acting on the light quark fields. They cannot easily be written 
in terms of an arbitrary T. For T = 7^, it is convenient to choose 

Q' 2 = e^' m q q^h Q , ,_ _ 

Q 3 = e^m q qV*Q, , # _ . <- (29) 

n , i<f/ _ <- Q 6 = e^ q(-w D)v^h Q . 

Q 4 = e 1<p q{-iD^)h Q , 

We neglect again operators whose matrix elements vanish by the equations of 
motion. A set of axial vector operators with the same coefficients bi(fj) is obtained 
from (28) and (29) by replacing q — > — q~7 5 and m q — > — m q . From tree level 
matching one computes 60 — b\ — 1 and b{ = for i > 2, while at leading 
logarithmic order (this may be extracted from ref. [27]) 

&oM = hQi) = z~ &IP , 

b 2 (f,) = -f -^/P + ^/P+l z -6/P f 

bM = f-fzW + lzW + l z -W, 

bM = 10 _ l z VP _ 2z -6/P ^ 

65^) = -f - + ^z-W - fz~W In 

W = "¥ + §zW-§z-W, 
where z = a s (mg) / a s (fj,) and (3 = 33 — 2Nf. 



(30) 



z 
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Let us now consider matrix elements of the effective current between a heavy 
meson and the vacuum. The application of the trace formalism to this particular 
case allows us to write 

(0\qTh Q \M(v)) = i^Tr[TM] , 

(0\qTiV^h Q \M(v)) = ^Tr[(FMv fl + F 2 ( f i) lfl )TM] , (31) 

( | q (-i^)niQ | M(v) ) = i Tr [ (F 3 (/x)^ + F 4 (/x) 7m ) T M] . 

Using an integration by parts and the equations of motion of the effective theory, 
one can show that [20] 

FiM = F 2 (/i) = FM = -^M [A _ Sm - m q ] . 

(32) 

F 3 (//) = -^[4(A-<H-m q ] • 

Additional parameters G\(fj,) and (^(/-O are induced by insertions of subleading 
operators from the effective lagrangian (7), in analogy to the functions 

It is clear from (32) that physical matrix elements depend only on the com- 
bination (A — 5m) and the form factors F,G\, and Gi, which must therefore be 
invariant under redefinitions of itlq. As an example, we present the result for the 
ratio of the decay constants of a heavy vector and pseudoscalar meson [20] 

V) [2&2M + &s(/i) + & 4 (/i) + W] - 16G 2 (A i )} . 

(33) 

From the fact that this ratio is /i-independent we conclude that the scale depen- 
dence of G2 involves m q and (A — 5m), in analogy to (26). 

The relations (32) allow us to give an operational definition of the combination 
(A — 5m), which is invariant under redefinitions of mg, without reference to the 
residual mass term. The reason is that this combination is induced when the usual 
covariant derivative acts on the light quark. In particular, it follows that 

(0| ^7 F ^' M ^^A-^. (34) 

Note that this relation expresses (A — 5m) directly in terms of properties of the 
light degrees of freedom, instead of in terms of a difference of heavy masses. We 
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may even propose (34) as a definition of the energy of the light degrees of freedom 
in the background of a static color source. This equation might also be useful in 
extracting the value of (A — 5m) from lattice gauge theory. 

5. 5. Summary 

We have demonstrated in detail that there is no important ambiguity in the 
choice of the HQET expansion parameter 1/mg, in the sense that all physical 
quantities are independent of this choice. This should hardly be surprising. The 
ambiguity is resolved by the introduction of a new quantity 5m, which is a non- 
trivial dynamical parameter of the effective theory. In order to be insensitive to the 
choice of the expansion parameter mn, hadronic form factors in the HQET must 
be defined in terms of matrix elements containing the operator YD^ = iD^ — Smv^. 
Physical quantities will then depend on these invariant form factors as well as on 
the combination (A — 5m), which replaces A in previous analyses. This combination 
is invariant under redefinitions of mq. 

While calculable in continuum perturbation theory, the residual mass term 
contains power divergences when the theory is regulated with a dimensionful cutoff, 
divergences which may be associated with incalculable nonperturbative processes 
[26]. The issue of whether such effects make the HQET useless beyond leading 
order in 1/mg is beyond the scope of this letter. It is important, however, to 
disentangle this larger question from that of whether the theory is plagued by a 
fundamental ambiguity in the choice of expansion parameter. We have shown that 
it is not, in the presence of the residual mass 5m. We believe that this is the 
correct framework within which to investigate further the possible importance of 
nonperturbative effects. 
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